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RAINDROP-SIZE DISTRIBUTIONS IN A HILLY REGION 
By C. A. NICHOLASS and P. R. LARKE 


SUMMARY 


Raindrop-size distributions have been measured in a hilly region of North Wales for the 
summers of 1972-74 and used to calculate empirical relationships between certain parameters 
important in radar and radio meteorology. 

The observed rainfall rates were between 0-1 mm h~? and 47-9 mm h-}, and about 34 000 
raindrop-size distributions were used to calculate various relationships by a linear regression 
technique. 

The drop-size distributions at the lower rainfall rates showed close agreement with the 
Marshall—Palmer exponential distribution, but for higher rates there was an increasingly 
noticeable deficiency of both very small and very large drops. 

The radar reflectivity factor (Z) and rate of rainfall (R) were calculated from each distribu- 
tion, and the empirical relationship of the form Z = AR, where A and B are ‘constants’. 
The constants calculated were higher than those measured by most other authors who used 
a similar regression technique. A comparison of areal radar measurements using the widely 
accepted value for the exponent (constant B) of 1-6 and the value of 1-8 derived from distro- 
meter data from North Wales produced an improvement in the accuracy of measurements. 
of about 10 per cent for rainfall rates less than 1 mm h-}, and 2-3 per cent for rainfall rates 
equal to or greater than 1 mm h-?. 

The Z-R relationships have been calculated for various synoptic types and storm by storm. 
These results show wide variability and it is concluded that radar systems should be calibrated 
using a few telemetering gauges in the area of interest. 

Empirical relationships between rate of rainfall, attenuation, median drop-size diameter 
and liquid water content have also been derived, and these confirm the results of other authors. 


INTRODUCTION 


Size distributions of raindrops have been measured and interpreted in many 
parts of the world, largely because they are important in determining the 
accuracy of measurements of precipitation using radar and assessing the 
problem of attenuation in radio communication caused by hydrometeors. 

It has been shown that the drop-size distributions are highly variable in space 
and time, and because of this the empirical relationships between the various 
parameters that can be derived from them (such as radar reflectivity factor, 
rate of rainfall and liquid water content) are subject to wide variations. (For a 
survey of these relationships see for example Stout and Mueller (1968).) 
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As a part of the Dee Weather Radar Project (Harrold et alii (1974)), raindrop 
spectra have been measured at two sites in North Wales to investigate their 
variability in a hilly region where accurate measurements of areal precipitation 
are difficult to achieve. Harrold et alii showed that although radar measure- 
ments of precipitation could be made with sufficient accuracy for many hydro- 
logical purposes, one of the major limitations to the accuracy obtainable was 
the change in the drop-size distribution over the area of measurement. The 
purpose of this paper is to describe the results derived from the very large 
number of drop-size spectra measured, to show the difficulties in achieving 
accurate rainfall measurements from a radar system, and to calculate the 
attenuation of microwave radiation caused by rainfall. 


INSTRUMENTATION 


The raindrop spectra were measured using a Joss distrometer* (Joss and 
Waldvogel (1967)). This instrument has a collecting area of 50 cm? and a crystal 
clock provides time pulses at one-minute intervals. The mechanical momentum 
of each raindrop is transformed into an electric pulse, and thereby the size of 
each drop is measured using a pulse-height analyser. Each drop is classified 
into one of 20 size categories ranging from 0-3 mm to 6-0 mm in diameter. 
However, it was found that in the first eight size intervals the numbers of drops 
measured alternated between many and a few for consecutive sizes, a defect 
attributed to slight inaccuracy in calibration. This was overcome by pairing 
size intervals, therefore reducing the first eight to four larger size ranges. 

The two sites chosen were at Carrog and Bala (see Figure I). The instrument 
was at Carrog during the summers of 1972 and 1974, and at Bala in the summer 
of 1973. The distrometer was mounted with its rim flush with the surrounding 
ground and close to a Plessey MM37 tipping-bucket rain-gauge which served 
as a check by providing an independent measure of rainfall. A comparison of 
gauge- and distrometer-measured rainfall is shown in Figure 2; the annual 
totals for the coincident periods are shown in Table I. 


TABLE I—RAINFALL MEASURED BY GAUGE AND DISTROMETER 


Year Gauge-measured Distrometer-measured Hours 
rainfall rainfall of data 
mm mm 
1972 104°0 IIg'l 1580 
1973 145°6 154°6 2256 
1974 93°6 82-2 1275 


Aniol (1974) indicated that a distrometer underestimates the rainfall in 
strong winds, an occurrence which would be expected in North Wales. Table I 
does not substantiate this but both Figure 2 and Table I demonstrate that, 
though errors are likely both with gauge and distrometer, the measurements 
are sufficiently alike for us to assume that the distrometer correctly measures 
raindrop distributions. 





* Dr Joss informs us that the word distrometer is derived from distribution meter and is the 
correct spelling for the instrument developed by him and Dr Waldvogel. The spelling 
disdrometer had previously been used for another instrument (Editor). 
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DROP-SIZE DISTRIBUTIONS 


The drop-size distribution can be defined by the summation over all sizes of 
the number of drops in a unit volume (Np) with diameters in the range D to 
D+ AD. 

Although the distrometer measures the number of drops (N) of given sizes 
hitting a sensor of area 50 cm? per minute (i.e. a flux) it is normal to express 
the drop-size distribution in terms of Np which is the number of drops per 
unit volume per unit size interval with a diameter D. The size intervals used 
and their respective raindrop terminal velocities are shown in Table II. 


TABLE II—TERMINAL VELOCITY OF RAINDROPS 


Category Minimum diameter Maximum diameter Terminal speed 
mm m s~* 
05 1-62 
o’'7 2°47 
3°36 
4°61 
5°41 
5°87 
6°40 
6:98 
7°48 
7:88 
8-21 
8-50 
8°77 
8-95 
9°05 
9°12 


© COI AUAWN 


66UNRK10 63 


CUHK 6IR HK OAK GUY 
ANAAWWWN NN Ee 


I° 
I° 
2° 
2° 
2° 
> 
3° 
3° 
4 
4 
5 


Marshall and Palmer (1948) related their drop-size distributions measured in 
Canada by Np=Noe—*” where No is the value of Np for D=o and A is a 
function of the rate of rainfall, A=41R-°?1 cm—!. The results obtained in 
North Wales are compared with this widely accepted distribution. 

Thecaiculated parameters have been compared using theempirical relationship 
in the form of the afore-mentioned Z-R relationship i.e. Z=AR® where A and B 
are ‘constants’. 

Stout and Mueller (1968), in their review, have shown that these ‘constants’ 
have in fact widely varying values dependent on location, synoptic weather type 
and probably the method of measurement. 

The mean raindrop-size distributions for a number of selected rates of 
rainfall have been computed for each of the three summers, by averaging the 
drop spectra for all minutes when the rainfall rate was within specific bounds. 
For example, all minutes with rainfall rates between 0-6 and 1-5 mm h™~? have 
been averaged to give the mean distribution for about 1 mm h7! (0-95 mm h7? 
in fact). If only minutes with exactly 1 mm h~? had been chosen, the data 
sample would have been rather small. The mean distributions for three fixed 
rates of rainfall are shown in Figure 3. For comparison the Marshall—Palmer 
(1948) distributions for the equivalent rainfall rates are also included. 

It can be seen that the mean distributions for each year were very similar to 
one another and to the appropriate Marshall—Palmer distribution. However, 
there was a tendency for fewer small and fewer very large drops to be measured 
than would have been expected from the Marshall—Palmer distribution. These 





TIVANIVY 
4O SALVU GALOIS ATWHL AOA SNOILNAIAISIG AZIS dowd NVAW AHI—€ gUnold 


(ww) g sajaworp dog 
02 St 
T T 








Meteorological Magazine, 105, 1976 








JOA@jU! @ZIS WW |. Jed @UNIJOA j1UN Jed Sdoup jo sequin 








JO [JOjuIDs UD@W ».y WW O69 9404 [jOjuIOI UDSW \-y WW 960 @jOJ |jOjulos UDaW )Yy WW 960 

QjO4 |JOyuIOs yO @Buos ».y WW S| ~-9-0 904 |jOyuros yo aBuos ;_y WW ¢-|-90 9404 |JOyuLos yo @BuDd4 ;-y WW ¢1-9.0 
DIOP sajnuiw 6907 DJOP sajnuiw {7g8E DIOP sajnuius $giZ 

¥Z6l DOWVD — — — — €Z6l VIV8 = —— — 2Z6l DOWVI = — — — 


904 |JOjUIOJ »_Y WW (| JaWjOYg PUD |}OYsIDW 

















1-Y WU 96°Z peor P[Noys oyeI [[eyUIes UvOUI ZL6] BOIIeD 


panuijuoo—€ aAnold 


(ww) q Jajawoip doig 


S@ 0% or" 
T | 














JOAsauU! @2IS WW [9 48d BWNjOA UN Jad sdoip jo saquiny 











GjO4 jjOyulos UDBW YY WW ZOE BJO |JOjUIO4 UdDSW |-Y WW ODE GjOJ [JOjulOs UOSW |-Y WW OO-E 

9404 OyuiOs yo Bus ».y WwW ¢-E-9.7 404 |jOjulos yO @Buos ;.Yy WU CE-97 JO |jOjulOs yO BBuos ;-y ww ¢-£-97 
OOP sajnuiw CE DIOP sajnuiw yy DIOP sajnuiw ~SE 

¥Z6l DORVI— — — — €Z46l VIVa-——— 7461 DOWVI—— — — 


3404 [poyuros t Y ww OF 4awjOg puo }o4ysow 














Meteorological Magazine, 105, 1976 





Meteorological Magazine, 105, 1976 


panuijuoa—t FANDLA 


(ww) g s2ewo!p dog 


of SZ 0z 
T 


St 
T 











————— aj03 jjoyuio) Uoaw , 
BOs YOyuIDs yo e6uo) 4 ww os-or 
DIOP Sonu 97 


41 SOWVI——-—- — 

















ww 705 —— 


aol yojuros uOew t y wu Lor 
9404 OyuIOs yO @Buos |.y WW O-¢~oy 
DIOP sajnuiw |Z) 


€i6l ViVE=-——— 





1DJ )Y WW OG JUDY PUo |jOYsSIOW 





uw SOS 


@j04 |jOjurOs yo @Buol ;-y WW O>-9P 


tZ6t 


saynuiw ~g 


oCawjes= 














JOAsaju) @z1S ww 1.0 49d awnjOA yun sad sdozp yo saquinyy 





Meteorological Magazine, 105, 1976 369 


measured deficiences become more marked at higher rainfall rates, as demon- 
strated in Figure 4, which shows the mean raindrop spectra for a number of 
rates from 0-95 mm h~? to 38-3 mm h~? at Carrog in 1972. (At the higher rates 
shown in Figure 4, data covering a very few minutes were included in the mean 
distributions.) The Marshall—Palmer distributions for 10 mm h~! and 40 mm 
h—? are also shown. 

It is worth noting that the maximum numbers of drops at higher rainfall 
rates occur at about 1-2 mm diameter. Peaks in the drop-size distributions 
have been observed by Caton (1966) using Doppler radar in England as the 
measuring instrument and by Hardy (1963) in the United States using a photo- 
electric spectrometer, although their maxima were for drop diameters between 
0-5 and 1-0 mm. 

Some authors (e.g. Joss and Waldvogel (1970)) have subdivided their data 
into synoptic types and found different distributions for each. The data from 
North Wales have been treated in this manner, by selecting the periods of 
rainfall ahead of warm fronts, cold fronts and occlusions, rain associated with 
depressions (called cyclonic) and shower rain (including thunderstorms) and 
then finding the mean drop-size distributions for given rainfall rates. It was 
found that for given rain types and for fixed rates the mean drop spectra were 
very similar in the three years. An example is shown in Figure 5 where the drop- 
size distributions for rainfall rates of about 3 mm h~! ahead of cold fronts are 
plotted for the three years, together with the Marshall—Palmer distribution for 
3 mm h-?. Once again the shortage of very small and very large drops can be 
seen. 

Figure 6 shows the mean distributions for four weather types (occlusions 
did not provide sufficient data) at about 1 mm h~! for Carrog 1972. It can be 
seen that for drops larger than 2 mm diameter, the four distributions are very 
similar. Below 2 mm there is a considerable difference between the types, in 
particular between frontal rain and showers. Surprisingly showers have more 
than twice as many small drops as warm frontal rain but only about 65 per cent 
as many medium-size drops. 


THE RELATIONSHIP BETWEEN RAINFALL INTENSITY AND RADAR REFLECTIVITY 
The parameters derived most frequently from drop-size distributions are the 
rate of rainfall R (in mm h~) and the radar reflectivity factor Z (in mm®m-*), 
It is widely accepted that an empirical relationship between them can be used 
to estimate rates of precipitation from radar measurements. However, the 
results of other authors around the world show that the relationship is highly 
variable and the use of a fixed equation could lead to serious error in the radar 
estimate of rainfall intensity. Np and D are known for each drop-size distribu- 
tion and since 


16 16 
Zo > Np,D® and Ro > Np,D?, 
i=1 i=1 


values of Z and R maybe calculated for each distribution. 

The relationship between Z and R is usually expressed in the form Z = AR®. 
However, since R is here regarded as the dependent variable, the least-squares 
fit linear regression is of log R on log Z. 

This procedure was used in the results which follow, the drop-size distribu- 
tions being measured by the distrometer over I-minute periods. Table III shows 
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the Z-R relationships in the three separate years calculated from log R on 
log Z linear regression for each minute when a rainfall rate of o-I mm h™ or 
more was measured. The Table also shows the number of minutes of data 
(that is, the number of drop-size spectra used in each year’s regression), the 
correlation coefficient, the range of measured rates of rainfall and the standard 
deviation of the constants A and B. 

Z-R relationships for different parts of the world have been quoted by 
several authors, but although a least-squares method was used it is not always 
stated whether a regression of log Z on log R or vice versa was calculated. 
Table III includes results by both methods for comparison. 


TABLE III—RADAR REFLECTIVITY (Z) AND RAINFALL RATE (R) RELATIONSHIPS 


Carrog 1972 Bala 1973 Carrog 1974 
log R log Z log R log Z log R log Z 
on on on on on on 
log Z log R log Z log R log Z log R 
Number of raindrop-size 
distributions 9814 15 764 8473 
Range of rainfall rates O'I to 47°9 Ol to 35°2 Ol to 34°4 
Correlation coefficient 0-940 0°936 0934 
A 


249 206 206 197 
Standard deviation of A 23 16 14 26 18 


B 1°74 1°53 85 1:60 1-76 1°52 
Standard deviation of B 0'026 0021 S 0:058 0:075 0:020 


It can be seen from Table III that the Z—R relationships for the three years 
are very similar, approximating to 250 R’® for a log R on log Z regression, 
and 200 R?® for a log Z on log R regression. (It was noted that about 1000 
minutes were needed to obtain a consistent Z—R relationship in each year.) 

There is considerable doubt which form of regression some past authors have 
used. The values of A and B in the Z = AR? relationship for North Wales 
calculated by a log Z on log R regression are very similar to those calculated 
by Marshall and Palmer (i.e. Z = 220 R*). However, calculating A and B 
from a log R on log Z regression results in higher values than those obtained 
by* most other authors who used this form of regression, though theirs are 
based on a relatively small number of drop-size spectra when compared with 
the very large number of distributions shown in Table III. 

It is apparent that, in order to improve the accuracy of measuring rainfall 
by radar methods over North Wales, a value of 1-8 for the exponent in the 
Z-R relationship should replace the value of 1-6 which is used at present (see 
Harrold et alii.) 

A comparison of results has shown that using a value of 1-8 the improvement 
in accuracy at low rainfall rates (<1-0 mm h7*) is about Io per cent. At rainfall 
rates of most interest (i.e. >1 mm h-) only a small improvement of 2-3 per 
cent is achieved. 

The lack of marked improvement is attributed to the fact that the radar 
samples aloft (about 300 m above ground at Carrog and 500 m above ground at 
Bala) whereas the relationships calculated in this paper are derived from drop- 
size distributions measured at ground level, which may have been affected by 
changes during the fall of the drops caused by low-level growth or evaporation, 
and size sorting by wind shears. 

It has been proposed (Joss et alii (1970)) that more accurate measurements 
of rainfall can be made using radar if a different Z—R relationship is used for 
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various synoptic weather types. Table IV shows the average Z-R relationships 
in the three summers for five synoptic weather types, together with the number 
of drop-size distributions used and the range of rainfall rates. 


TABLE IV—RADAR REFLECTIVITY (Z) AND RAINFALL RATE (R) RELATIONSHIPS 
(FROM LOG R ON LOG Z REGRESSION) FOR VARIOUS SYNOPTIC TYPES(Z=AR?) 


Carrog 1974 
A B 


Carrog 1972 Bala 1973 
A B A B 


Cold fronts 224 1°72 268 1°84 220 1-90 


Number of distributions 
Maximum rainfall rate* 
Warm fronts 

Number of distributions 
Maximum rainfall rate 
Occlusions 

Number of distributions 


Maximum rainfall rate 
Cyclonic rains 

Number of distributions 
Maximum rainfall rate 
Showers 

Number of distributions 
Maximum rainfall rate 


264 


422 


211 
(only one storm) 


2477 
12°3 


778 
123 


2°0 


479 


1:60 


201 


3814 
74 
229 1°79 


25°4 
194 
805 


12°8 

232 1:70 
2613 
24°9 

254 1-79 
942 
35°2 


260 


316 


1736 
III 


I° 
1886 


21°4 


I° 
1716 


34°4 


* The minimum rainfall rate is 0-1 mm h-}. 


The results show that by reducing the data samples by subdivision there is 
a greater uncertainty in the Z—R relationship, and that any increase in accuracy 
which might be achieved would not justify the extra effort involved. 

An attempt was made to relate the ‘constants’ A and B for the different 
synoptic types to low- and medium-level wind velocity but little correlation was 
established. 

Table V shows the storm-to-storm variation of Z-R relationships, using 
Carrog 1972 as an example. A few storms with very small ranges of rainfall 
rates have been excluded. 

It can be seen that if only a few storms, or parts of storms, had been used as 
the data sample on which to base the Z—R relationships a large error would 
be possible. Furthermore, if the mean Z-—R relationship were to be applied to 
radar rainfall measurements, considerable day-to-day inaccuracy would result 
unless the radar system was additionally calibrated using a few widely spaced 
telemetering rain-gauges in the areas of interest (see Harrold et alii). The use 
of a number of distrometers to provide real-time Z—R relationships as an input 
to the radar would not seem to be feasible, because of the large amount of data 
processing involved. 

Mean relationships between Z and R for all synoptic types for the three 
summers are illustrated in Figure 7. 


THE RELATIONSHIP OF ATTENUATION TO RATE OF RAINFALL 


Many authors (see for example Harrold (1974)) have stressed the importance of 
attenuation, caused by precipitation particles, on the transmission of radio 
signals at millimetric wavelengths. The ability of radar to measure precipitation 
is affected to a lesser extent by attenuation but nevertheless must be considered. 

Gunn and East (1954) comprehensively reviewed the physics of attenuation 
and showed that it was derived from a complicated function dependent on the 
wavelength of the radiation, the temperature and the drop-size distribution. 
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TABLE V—STORM RADAR REFLECTIVITY (Z) AND RAINFALL RATE (R) 
RELATIONSHIPS AT CARROG 1972 


Date Number of Z-R relation Correlation Range of rain- 
drop-size coefficient fall rates 
distributions B (mm h-?) 
Cold fronts 
29 Apr. 262 0925 orI- 
25 May 62 0-964 
2 June 287 1°85 0°949 
8 June 1°51 0-976 
27 June 1-91 0940 
4 Aug. 1-92 0°944 
5 Aug. 1°65 0°977 
8 Aug. I-31 0°980 
17 Aug. 1°58 0°959 
6 Sept. 1:78 0°945 
8 Sept. 1-61 0'906 
11 Oct. 1-64 0'944 
Warm fronts 
28 Apr. 0°952 
2 June 1-32 0°975 
7 Aug. 1°43 0-962 
8 Sept. 0°967 
Occlusions 
20 June : 0°906 
Cyclonic rains 
26 May 1°65 0°959 
9 June 255 320 1°52 0-961 
1 Aug. 347 173 I°51 0°970 
8 Aug. 368 119 1°51 0-958 
Showers and 
thunderstorms 
May 407 232 1°65 0°953 
1 June 139 224 1°55 0°967 
27 June 223 187 1°83 0917 
6 July 67 293 1-76 0°903 
29 July 29 349 1°59 0°985 O:I-16°7 
30 July 61 487 1°51 0961 OI— 3°5 
9 Aug. 45 241 1°71 0-948 OrI-16°1 
13 Aug. 68 165 1°59 0-970 OI-— 33 
7 Sept. gI 204 1-72 0'924 OI— 38 
13 Sept. 164 237 1-7I 0953 O:I-I1°7 


They showed that at wavelengths which are long compared with the drop 
diameter the one-way attenuation (K), (i.e. half the total attenuation caused 
by precipitation en route to target and on return), in decibels per kilometre is 
given by 


K = 0°4343 2Qt, 

where the summation is over one cubic metre and Q; (the total absorption 
coefficient) is in square centimetres. Herman et alii (1961) have produced tables 
of normalized attenuation at 0o°C (=4Q;/7D7*) for a variety of wavelengths 
and drop sizes. From this quantity it is possible to calculate the attenuation 
caused by the measured drop-size distribution and hence to relate this to the 
rate of rainfall or radar reflectivity, both of which are derived from the distribu- 
tion. 

For each of six wavelengths frequently used in radar or radio work, the 
attenuation has been calculated in each minute when the rate of rainfall was 
o°I mm h~ or more. A mean relationship of the type 

K = ER’ 
has been found by a least-squares linear regression of log K on log R. Table VI 
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shows the results of these regression calculations and a comparison with those 


calculated by Wexler and Atlas (1963). In every case the correlation coefficient 
was 0:97 or higher. 


TABLE VI—ATTENUATION (K) AND RAINFALL RATE (R) RELATIONSHIPS BASED 
ON ATTENUATION CONSTANTS E,F AT 0°C 


Wavelength Carrog 1972 Bala 1973 Carrog 1974 Marshall—Palmer 
cm after 
Wexler-Atlas 
E F E F E F E 


10°0 0:00071 0:98 0:00075 0°95 0°00070 0°99 0:00082 
55 0:0026 0:94 0:0028 0°89 0:0026 0-92 0°003I 
3°21 O-oll 1-Ol O-OIl 0:97 O-Ooll 0-98 0-013 
1:87 0:038 I'll 0:041 I-10 0°037 I-12 0:050 
1:24 0:099 1:05 O'105 1:03 0°097 1:05 013 
0°86 O'215 1:07 0:229 1:05 0-214 1:06 0°31 
The attenuation constants were computed by Herman et alii using the Marshall and Palmer 
drop-size spectrum. 


It can be seen that in each year the two ‘constants’ in the K—R relationships are 
very similar. At the wavelengths usually used for radar measurement of 
precipitation (5:5 to 10 cm) the attenuation is very small even at high rainfall 
rates. Using a 5:5-cm radar the attenuation at rainfall rates of 1 mm h™ 
and 5 mm h~ which are typical of North Wales result in attenuation of 00026 
dB km“ and 0-0118 dB km“. 

The relationships derived from each year are also similar to those derived 
using a Marshall—Palmer distribution. The differences are attributed to the 
mean drop-size distributions described earlier showing slight variations from 
the Marshall—Palmer distribution at most rainfall rates. 

The K-R relationships have been calculated using the full range of rainfall 
rates, although a high proportion of the data are for relatively low rates. For 
example, at Carrog in 1972, out of the 9814 drop spectra used, 9565 were for 
rates of rainfall of 5-o mm h~ or less. Thus the mean attenuation relationships 
obtained are most applicable at low rainfall rates, although the attenuation 
is more important at the higher rates, which fortunately occur infrequently. 
In an attempt to show how the attenuation relation is dependent on rainfall 
rate, the data have been divided into three overlapping ranges of rainfall rates, 
namely 0-I-—5°0, I-O-10°0 and greater than 5:0 mm h~’, and the K-R relation- 
ships calculated by linear regression for drop-size distributions within these 
ranges. The results in the three years are very similar again, so only those for 
Carrog in 1972 are shown in Table VII. 

Table VII shows that at wavelengths of 10 cm and 0-86 cm there is little to 
choose between the three regression equations. For the other wavelengths it can 
be seen that there is a tendency for the regression equations to have steeper 
slopes at higher rainfall rates, implying that the attenuation may become a 
more serious problem at very high rates than would be suggested by the original 
relationships shown in Table VI. A very large number of drop-size distributions 
are required, particularly at higher rates before the attenuation relationship 
can be defined with confidence at the rates of rainfall when it is more important. 

Empirical relationships between attentuation (K) and radar reflectivity (Z) 

can be derived from the previously calculated Z—-R and K-R relationships. 
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For example, at Carrog in 1972 for a wavelength of 5-5 cm 
Z = 249R'4 and K = 0:0026R°™ 


lead to the relationship 
K = 1°32 x 107* 27°. 


TABLE VII—ATTENUATION (K) AND RAINFALL RATE (R) RELATIONSHIPS OF 
FORM K=ERF as A FUNCTION OF R AT CARROG 1972 


Wavelength Range of Number of K-R relationship Correlation 
rainfall rates distributions coefficient 
cm mm h-* E F 
10°0 OI- 5:0 9565 0:00071 0-98 0-978 
1:0-10°0 2575 0°00072 0:90 0'971 
More than 5:0 249 0:00067 0:95 0-992 
55 O'I- 5:0 0:0026 0°93 0°973 
1:0-10°0 As above 0:0025 0°95 0-976 
More than 5:0 0:0020 I-10 0:973 
3°21 OI- 5:0 0-0106 1-00 0-980 
10-10°0 As above 0-0102 1:06 0°974 
More than 5:0 0:0074 1-23 0956 
1:87 OI- 5:0 0:0383 I-ll 0°985 
1:0-10'0 As above 0°0380 1-09 0°978 
More than 5:0 0:0317 1-18 0974 
1°24 O'I- 5:0 0:0983 1:05 0986 
I'0-10'0 As above 0:0966 1:07 0°983 
More than 5:0 00815 1-16 0-980 
O'I— 5:0 0°216 1-07 0'987 
10-100 As above 0°215 1-06 0-986 
More than 5:0 0*209 1-07 0°990 


OTHER METEOROLOGICAL PARAMETERS DERIVED FROM DROP-SIZE DISTRIBUTIONS 


Although they are not frequently used in connection with radar rainfall measure- 
ments two other parameters, the liquid water content (M) and the median- 
volume drop diameter (Do) can be calculated from the drop-size distribution. 

The liquid water content is the mass of water (in grams) contained in 1 m$ 
and is calculated by the summation of the masses of water contained in each 
drop category in 1 m°; 


that is M = > NpMp AD, 
1m* 


where Mp is the mass of a drop with a diameter in the range of Dto D+ AD. 

The liquid water content is related to the rate of rainfall by an empirical 
relationship of the type M = GR”. The ‘constants’ can be derived by a least- 
squares linear regression of log M on log R, as in previous sections. 

The liquid water content has been calculated for each drop-size distribution 
in the three years and the results of the linear regressions are shown in Table VIII. 
Also shown are the ‘constants’ in the relationship Z = JM? which can be 
derived from the Z—R and M-R relationships. 

The median-volume drop diameter (Do) has been defined by Atlas (1953) as 
that drop diameter which divides the drop distribution in such a way that half 
the liquid water content is contained in drops greater than Do. As the liquid 
water content is related to the rate of rainfall by an empirical power relationship 
as above, it is reasonable to seek a similar relationship between the median- 
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TABLE VIII—-LIQUID WATER CONTENT (M/), RAINFALL RATE (R) AND RADAR 
REFLECTIVITY (Z) RELATIONSHIPS 


Site Number of M=GR* relationship Correlation 2Z=JM? relationship 
distributions (log M on log R) coefficient (direct calculation) 
G H J P 


Carrog 1972 9 814 0:0794 0°93 0-960 2°79xX 104 1°86 
Bala 1973 15 764 0:0863 0:90 0946 3°93 X 104 2°05 
Carrog 1974 8 473 0:0805 0°93 0:958 2°87xX 104 1:90 


volume drop diameter and rainfall rate. The results of log Do on log R regres- 
sions for the three years’ data are shown in Table IX. Atlas reported that the 
scatter about a line of the form Do = SR’ (where Dp is in mm and R in 
mm h~?) is fairly large, and this is demonstrated by the rather low correlation 
coefficients in the Table. 


TABLE IX—MEDIAN-VOLUME DROP DIAMETER (Do) AND RAINFALL RATE (R) 
RELATIONSHIPS 


Site Number of Do=SR’ relationship Correlation 
distributions Ss ‘i coefficient 
Carrog 1972 9 814 0°994 0°20 0-630 
Bala 1973 15 764 0-961 0°24 0-677 
Carrog 1974 8 473 0-962 o'21 ; 0°623 
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THE HEIGHT OF THE FREEZING LEVEL DURING RAINFALL 


OVER THE BRITISH ISLES 
By C. G. COLLIER 


SUMMARY 


An analysis of the height of the freezing level over the British Isles is presented for November, 
January and March using the ten years of data 1964-73. The relationship of the spatial and 
temporal variations of this parameter to twelve-hour rainfall amount is investigated. 


I. INTRODUCTION 


Recent work on the measurement of rainfall using radar (see for example 
Harrold et alii (1974)) has demonstrated that radar rainfall measurements may 
be made with an accuracy of 15 to 20 per cent over periods of three hours, 
provided that the melting level is not intersected by the radar beam. When the 
melting level intersects the radar beam the accuracy is significantly decreased, 
since the presence of melting snow in the beam causes the power received at the 
aerial to be greater than that from rain at the same rate of precipitation (see 
Battan (1973)). Methods which attempt to correct the radar rainfall estimates 
when this occurs have been developed (see Harrold and Kitchingman (1975)). 
Nevertheless the errors that the melting level may introduce into the radar 
estimates cannot be completely removed. In order to examine how often this 
error may affect radar measurements in any part of the British Isles, a study 
is presented of the spatial variation of the height of the freezing level meaned 
over a ten-year period. The study is only concerned with the o°C level. The 
centre of the melting layer is below this height, occasionally by several hundred 
metres (see Battan (1973)). 

Previous work on the height of the freezing level by Murray (1950) has 
indicated its seasonal variations over the British Isles using six years of data 
1942-47. The results are not related to rainfall, and therefore not directly 
applicable to the estimation of frequency of error occurrence due to melting- 
level effects in radar rainfall measurements. 

In processing the data, when the temperature at the surface was less than 
o°C the freezing level was taken as being at the surface. For cases in which a 
sounding had two or more heights at which the o°C isotherm was intersected, 
only the lower height was considered, as this height is closest to the true height 
of the melting level relevant to radar measurement accuracy, as explained 
above. The freezing-level heights have been extracted from the radiosonde data 
recorded twice daily (00 and 12 Gr) at each of the routinely manned British 
Isles upper-air stations. The rainfall data have been extracted from the Daily 
Weather Reports for surface stations within 25 km of the upper-air stations. 
Twelve-hour rainfall totals spanning the times of the radiosonde ascents (the 
‘rainday’ begins at 09 GMT) have been used in the analysis. Since the atmospheric 
structure varies considerably with time, the use of one ascent to represent 12 
hours of rain imposes a severe limitation on the analysis, but more frequent 
upper-air temperature data are not available as routine. 
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2. SPATIAL VARIATION OF FREEZING-LEVEL HEIGHT IN WINTER 


The variation of freezing-level height over the British Isles has been examined 
by Murray (1950) for each month of the year covering the period 1942-47. 
The ten years of data 1964-73 considered in the present study give mean 
freezing-level heights, use being made of all the data for the months of Novem- 
ber, January and March, which agree well with those given in the earlier study. 
During the winter months the mean height of the freezing level, all data being 
considered without regard to rainfall amount in 12 hours, decreases uniformly 
from approximately 1100 metres over the south-west of the British Isles to 
approximately 600 metres over the northern tip of the Shetlands, Figure 1. 
In summer this lowering of freezing-level height is from the south-east to the 
north-east (see Murray (1950)), and the transition from summer to winter is 
illustrated by the month of November shown in Figure 1. The standard devia- 
tions about the mean values are between +300 metres and +800 metres, see 
Figure 2. 

It is of interest to potential users of a radar system measuring rainfall to 
examine whether or not these mean maps are appropriate for all rainfall rates. 
The lower the height of the freezing level, the closer to the radar site its inter- 
section with the beam, and therefore the greater the area of radar coverage 
affected by measurement errors introduced by its presence. Several ranges of 
12-hour rainfall amount have been considered. 


2.1 Rainfall rate zero or 0°05 mm in 12 hours 

The distributions of mean freezing-level heights for zero or negligible rainfall 
in 12 hours are shown in Figure 3. Compared with the mean distribution for all 
rainfall rates (Figure 1), the freezing level is higher over the south-west of the 
British Isles in November, lower over the north-east of the British Isles in 
January, and slightly higher over the south-east of the British Isles in March. 
Table I gives the percentage frequency of occurrence of the freezing level below 
specified heights for three rainfall-rate ranges. This Table gives an indication 
of the amount of data available for each rainfall-rate classification. About 50 
per cent of the data lie within this rainfall-rate range (i.e. 0-0-05 mm/I2 hours). 


2.2 Rainfall rates between 0:05 mm in 12 hours and 5 mm in 12 hours 

In November and March the variation of mean freezing-level height for this 
class (Figure 4) resembles quite closely the mean pattern for all rainfall rates, 
but in January it is significantly different, showing a marked rise from the south. 
In January over Cornwall the mean height of the freezing level for this rainfall 
rate range is about 200 metres higher than given by the mean pattern (Figure 1), 
and this height increase is reflected throughout the British Isles. Approximately 
45 per cent of the data fall within this rainfall-rate range (Table I). 


2.3 Rainfall rates between 5 mm in 12 hours and 10 mm in 12 hours 

Figure 5 shows a lowering in the mean height of the freezing level for this 
rainfall-rate range in both November and March compared with the distribution 
for all rainfall rates (Figure 1), in November from the north-west, and in March 
from the north over the eastern part of the British Isles. This effect is most 
marked in March when it produces @ strong gradient of freezing-level height 
from east to west over England, and from north-east to south-west over Scot- 
land. In March the mean freezing-level heights for this rainfall-rate range are 
50 per cent lower than the mean values for all rainfall rates (Figure 1) over the 
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eastern part of the British Isles. In January the tendency for the freezing level 
to rise is continued from the south-west, producing a mean pattern which 
changes by 700 metres across central England. Table I indicates that rainfall 
rates in this range make up less than Io per cent of the data sample. 

Since only a very small percentage (less than 3 per cent) of cases have rainfall 
rates greater than 10 mm/I2 hours, it is felt that the data sample for these 
rainfall rates is too small to be discussed with any confidence. 


3. DISCUSSION 


The freezing-level mean height patterns show large differences from one rainfall- 
rate range to another, but the significance of the changes depends upon the 
frequency with which the rainfall rates occur. Figure 6 shows the percentage 
frequency of 12-hour rainfall for Ringway Airport, Manchester (taken with the 
upper-air data for Aughton, Liverpool) for the period 1964-73. Fewer than 
2 per cent of the data are composed of cases in which the 12-hour rainfall is 
greater than 10 mm, 4 per cent (of data) have rainfall in the range greater than 
5 mm to 10 mm, and 60 per cent (of data) have zero rain or a trace (Section 2.3). 
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FIGURE 6—PERCENTAGE FREQUENCY OF I2-HOUR RAINFALL IN JANUARY FOR 
RINGWAY AIRPORT (TAKEN WITH THE UPPER-AIR DATA FOR AUGHTON), 1964-73 


The marked differences from the mean patterns for the higher rainfall 
amounts may be related to the different frequencies of synoptic weather types 
through the winter months. However, there are only a small number of cases 
with rainfall rates in these ranges, and therefore the data sample is too small 
for any firm conclusions to be drawn from it. The frequency of daily weather 
types over the British Isles has been discussed by Lamb (1953, 1965), although 
the weather types are defined without reference to rainfall. Murray and Lewis 
(1966) extended this work, and suggested a system of indices which could be 
used to indicate the broad-scale circulation near the British Isles. In later papers 
Murray (1968, 1972) demonstrated the relationship of these indices to monthly 
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rainfall over England and Wales, and to monthly mean temperatures in central 
England. The tendency in winter months for low (high) rainfall to be associated 
with low (high) temperatures, and the tendency in the summer months for low 
(high) rainfall to be associated with high (low) temperatures were also discussed 
(Murray (1968)). This may account for the lowering of the freezing level in 
March (compared with the distribution for all rainfall rates (Figure 1) observed 
for the high rainfall-rate range (Figure 5) as March may be regarded as a transi- 
tion month from winter to summer. The freezing level rises (Compared with the 
distribution for all rainfall rates (Figure 1) for the high rainfall-rate ranges in 
January, a winter month. 

It is not the purpose of this paper to seek detailed correlations between 
rainfall, height of the freezing level and synoptic weather type as represented 
by the indices by Murray and Lewis (1966). To do so would require a much more 
detailed discussion than has been presented in this paper. The above represents 
a suggestion where the explanation of the form of Figures 3 to 5 may lie. 

Table I gives the percentage frequency with which the o°C level is below a 
specified height, for a particular rainfall rate. The effects of latitude and season 
are apparent in the Table. Figures 3 to 5 and Table I may be used to estimate 
the mean distance, for any radar site in the British Isles, at which the radar 
beam intersects the freezing level, for various ranges of 12-hour rainfall amount 
in winter. For a radar situated 100 m above sea level with a 1° wide beam and 
the beam axis elevated by 0°5° from the horizontal, the 0°C level intersects the 
beam at about 50 km range from the radar when it is at a height of 1 km. The 
stippling on Figures 3 and 4 shows the areas where the 0°C level would intersect 
the beam of this radar at ranges less than 50 km. 

Figure 7 indicates the mean percentage number of occasions with the 12-hour 
rainfall amount in the range 0:05 mm to § mm, when the 0°C level would inter- 
sect the beam of a radar, elevated by 0°5° from the horizontal, within 50 km 
of the radar site. The radar is assumed to be situated 100 m above sea level at 
each of the radiosonde sites, and at each site on the maps the upper figure is the 
average number of occurrences in the month with rainfall in the range 0-05 mm 
to 5 mm in 12 hours. The lower figure at the stations is the average number of 
occurrences in the month with rainfall in the range 0-05 mm to § mm in 12 hours 
when the o°C level intersects the radar beam within 50 km of the radar site. 
These maps enable the mean percentage amount of data, in the rainfall range 
0-05 mm to § mm in 12 hours, likely to be effected by measurement errors intro- 
duced when the o°C level intersects the radar beam within 50 km of the radar 
site, to be derived for any part of the British Isles. Ten years of data have been 
used to produce Figure 7, 1964-73, and it must be stressed that any one year 
may differ markedly from average. For example, for the C-band Llandegla 
radar used in the Dee Weather radar Project (see Collier et alii, 1975) (400 m 
above sea level) in January 1974, 61 per cent of the radar data collected were 
affected by measurement errors introduced by the presence of the melting level 
in the beam within 50 km of the radar site. However, in March 1974, only 
28 per cent of the radar data were affected. The corresponding mean figures 
derived from Figure 7 are 50 per cent in January and 60 per cent in March. 

Given the frequency distribution of the rainfall-rate ranges near a particular 
radar site, the mean temporal variations of the range at which the radar beam 
intersects the o°C level may be obtained. For example, for the Llandegla C-band 
radar, the data given in Table I for Aughton may be used assuming a beam 
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elevation of 0-5°. The mean freezing-level height over North Wales for rainfall 
rates between 0:05 and 5 mm ‘in 12 hours in November is 1250 metres, in 
January is 1100 metres, and in March is 1000 metres. On average, rainfall rates 
in this range occur on 7 days per month (corresponding approximately to a 
percentage probability of 35 per cent) in winter, and on 4 of those days or less 
the freezing level will be intersected at a distance of about 41 km (corresponding 
to a freezing-level height of 1250 metres which occurs with a percentage proba- 
bility of about 20 per cent for this rainfall-rate range) in November, 30 km in 
January and 27 km in March. However, for a rainfall-rate range of 5 mm in 
12 hours to 10 mm in 12 hours, the mean heights are 1300 metres in November, 
1100 metres in January and 950 metres in March. On average, rainfall rates 
within this range occur on only 1 day per month or less, and on half of this day 
the freezing level will be intersected at a distance of about 45 km in November, 
36 km in January and 30 km in March. 


4. CONCLUSION 
The variation of the height of the freezing level over the British Isles in winter, 
shown in Figure 1 for all rainfall rates, reflects the dependence of this parameter 
on latitude and season. However, there is a relationship between mean freezing- 
level height and 12-hour rainfall amount. 

The data presented in this paper allow mean estimates to be made of the 
distance from a radar site at which the freezing level intersects the radar beam 
in winter. The percentage amount of radar data with 12-hour rainfall in a given 
range affected by melting level errors may also be estimated. 
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